The Evolution of the Large-Scale Tail of Primordial Magnetic Fields 



Karstcn Jcdamzik 
Laboratoire de Physique Theorique et Astroparticules, UMR5207-CNRS, 
Universite Montpellier II, F-34095 Montpellier, France 



o 
o 

Q 

(N 

o 

u 

o 



> 



Giinter Sigl 

II. Institut fur Theoretische Physik, Universitdt Hamburg, 
Luruper Chaussee 149, D-22761 Hamburg, Germany 

Cosmic magnetic fields may be generated during early cosmic phase transition, such as the QCD- 
or electroweak- transitions. The magnitude of the remainder of such fields at the present epoch 
crucially depends on the exponent n of their (initially super-Hubble) large-scale tail, i.e. B\ ~ A~". 
It has been claimed that causality requires n = 5/2, contrary to much earlier claims of n = 3/2. 
Here we analyze this question in detail. First, we note that contrary to current belief, the large- 
scale magnetic field tail is not established at the phase transition itself, but rather continuoulsy 
evolves up to the present epoch. Neglecting turbulent flows we find n = 7/2, i.e. very strongly 
suppressed large-scale fields. However, in the inevitable presence of turbulent flows we flnd that 
the large-scale magnetic field tail has sufficient time to evolve to that of the fluid turbulence. For 
white noise fluid turbulence this yields n = 3/2 up to a certain scale and n = 5/2 beyond for the 
magnetic fleld spectrum. This picture is also not changed when primordial viscosity and fluid flow 
dissipation is taken into account. Appreciable primordial magnetic fields originating from cosmic 
phase transitions seem thus possible. 



I. INTRODUCTION 

The origin of galactic- and cluster- magnetic fields is 
still unknown (cf. Q). For long, one viable possibility 
was the generation of magnetic fields during some early 
magnetogenesis epoch which than could act as seed 
fields for a galactic dynamo or even under optimistic 
conditions provide the complete magnetic field after grav- 
itational collapse Q. Another viable option is the gen- 
eration of magnetic fields in stars and the expulsion of 
these fields into the intergalactic medium after the end 
of the life time of stars Q . In either scenario many open 
questions remain. Irrespective of the origin of galactic 
magnetic fields it is interesting to know if substantial 
primordial relic magnetic fields could have survived the 
evolution of the early Universe. 

This is particularly so due to recent claims of a lower 
limit on the magnitude of magnetic field strength in the 
intergalactic medium oi 10^^^ — lO^^^Gauss (depen- 
ing on coherence length scale) Q. These claims are 
due to 7-ray observations of distant blazars by Fermi 
and HESS and detailed analysis of 7-ray showers on the 
extragalactic- and cosmic microwave- background radi- 
ation in the presence and absence of intergalactic mag- 
netic fields (cf. 0). They seem to indicate that a large 
volume fraction / ^ 0.5 of voids a filled by magnetic fields 
which seems only difficult to achieve in scenarios of out- 
flows/winds from galaxies. Q 

It has been claimed that phase transition gener- 
ated primordial magnetic [lOl fields are too weak to be at 
the origin of the observed galactic- and cluster- magnetic 
fields, even under very optimistic assumptions on galactic 
dynamo amplification. The only exception would be here 
the case of completely helical fields produced during the 
QCD transition, nevertheless, requiring substantial dy- 



namo amplification as well The main point here 

is that due to "causality" arguments, it is claimed, the 
spectral index n in 



1 

8^ 



d^fc'BtBk = e"^ 



2n 



(1) 



would be large n 5/2 yielding a very blue spectrum and 
not much power on scales much above the phase transi- 
tion Hubble radius. Here Bk is the Fourrier transform 
of magnetic field, k is wave vector, and is magnetic 
energy density. 



In this paper we analyze the question of the spectral 
index in detail. In the next subsection wc compute the 
magnetic field tail arising from the superposition of un- 
correlated mag netic dipoles in vaccum, exemplifying that 
n = 3/2 [l2| is possible even for supper- Hubble dis- 
tances. In Section HI we replace the vacuum by the 
highly conductive medium in the early Universe, find- 
ing a suppressed n = 7/2. In Section IV we include the 
realistic existence of fluid flows in our analysis analyzing 
the evolution of the large-scale tail first, in the turbu- 
lent fiuid, and then in a fluid described by alternating 
periods of turbulent and viscous magnetohydrodynamics 
(MHD), as shown to occur in Ref [13|. In both cases we 
find n = 3/2 up to some large scale. Conclusions are 
drawn in Section V. 
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II. THE MAGNETIC FIELD SPECTRUM DUE 
TO UNCORRELATED DIPOLES 

For a magnetic dipole density M, in the magneto-static 
approximation the magnetic field can be written as 

^ I ^3^,3(r-rOM(rO.(r-rO M(rO|r-rr 
J |r-r'|5 

where /xq is the magnetic susceptibility in vacuum. From 
this one can easily show that 



(BkB_k) cx / d\' (M(r)M(r')> e 



(2) 



For uncorrelated dipoles, (M(r)M(r')) cx 5{y — r'), this 
gives (BkB_k) =const and therefore, comparing with 
Eq. ([1]), a slope esik) cx k^, thus n = 3/2. The magneto- 
static approximation is good at scales much smaller than 
the light-travel distance and thus certainly at Mpc scales, 
which are much smaller than the Hubble radius today. 



III. MAGNETIC FIELD EVOLUTION 
NEGLECTING FLUID FLOWS 

We assume that the total electrical current density con- 
sists of an external source jox and an Ohmic current in- 
duced by the electro-magnetic field. 



(t(E + v X B) , 



(3) 



with a the conductivity. Here the non-vanishing exter- 
nal current jox is assumed to be dynamically forced dur- 
ing the epoch of magnetogenesis. Upon neglecting the 



displacement current 5tE for the total current one has 
Amperes law 47rj = V x B. Writing the external current 
as the rotation of the dipol-density M, 



jcx = V X M , 



(4) 



inserting E from Eq. Q into Faraday's induction equa- 
tion 



dfB = -V X E 



(5) 



and using V • B = one has 



9fB = V X (v X B) + ?7 [AB - 47r (AM - V (V • M))] , 

(6) 

where we have used the resistivity 77 = l/(47r(T). 

We now express these equations in the expanding Uni- 
verse in terms of redshift z and at the same time perform 
a Fourier transformation of the spatial dependence. Us- 
ing = + z)H{z)]dt with H = a/ a the Hubble 
constant and the scale factor a = l/(l + z) and neglecting 
the fluid velocity term in Eq. ([6|) one then obtains for the 
evolution of the Fourier component Bk(2) for comoving 
wave number k. 



a.Bk(z) 



l + z 



(7) 



k2 [Bk(z)-47rM^(z)] , 



where Mj|r(z) = Mk(2;) — (k-Mk(z))k is the component 
of Mk(2:) perpendicular to k, with k the unit vector in 
direction of k. The solution of Eq. ([7]) which vanishes 
before the source term is switched on is given by 



J 



B(k,z) = 47rk2(l 



cxp 



-k^ 



H{z) 



77(z')M_L(k, z') 
{l + z')H{z') 



dz' . 



(8) 



It is convenient to define a critical wave number as 

-1/2 



kr{z) 



H{~z) 



dz 



(9) 



For k kr the exponential in Eq. ([5]) cuts off contribu- 
tions from z + Az > z + H{z)/[k'^{l + z)j]{z)] so that one 
obtains 



Bk(z)~47rM^(z). 



(10) 



This is the small-scale instantaneous limit in which the 
field is given by the superposition of dipole fields which 
are only subject to redshift since the second term in 
Eq. ([7]) approximately vanishes. In the opposite limit. 



k <^ kr, the exponential in Eq. ^ can be set to one and 
one obtains 

Bk(z) ^ 47rk(l + z) ^-—^^^dz 



4n 



kr{z) 



M^(z), 



where in the second estimate we have used that 



M^(z') 



1 + z' 
l + z 



mi{z) 



(11) 



(12) 



In this limit the resulting field is thus frozen into the 
plasma and strongly suppressed compared to the source 
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term by a factor {k/krY due to the screening by the con- 
ductivity of the medium. It is easy to see that the critical 
wave number becomes smallest at late times, when elec- 
trical conductivity and Hubble constant become smallest. 
Using the Spitzer resistivity rj ~ Trmy^e^/T'^/^ one 
finds 



T 



1/2 



fcr(r) ~ 3 X 10^ j pc-ifor T < leV, (13) 

where Tq is the present day CMBR temperature. This 
should be compared to, for example, the comoving QCD 
Hubble radius wave vector, fcqcD ~ 6pc~^(T/100 MeV), 
showing that screening of the plasma is extremely effi- 
cient, leading only to highly suppressed magnetic fields 
on large scales. The supression here, n = 7/2 is sub- 
stantially larger that that claimed to apply in Ref. 
However, we will see in the next section that the con- 
clusions change drastically when fluid motions arc taken 
into account. 



IV. MAGNETIC FIELD EVOLUTION WITH 
FLUID FLOWS 

As post-inflationary "causal" magnctogcncsis in the 
early Universe most likely takes place during cosmic first- 
order phase transitions, and during such transitions the 
generation of cosmic turbulence is inevitable, the com- 
plete neglect of fluid motions is thus unrealistic. In 
the one-fluid approximation of magnetohydrodynamics 
(MHD) the velocity field v evolves according to 



^fv + (v • V)v = 



(V X B) X B 



(14) 



where p is the fluid density and fu is the viscous damping 
force [13 ■ We assume now that the fluid kinetic energy 
density (at least initially) dominates magnetic energy 
density. Furthermore we assume that the phase tran- 
sition induced magnetization M has led to the build-up 
of some magnetic seed field on relatively small scales (cf. 
Sec. HI), but nothing on large scales. We are then in- 
terested to know how much magnetic field will be gen- 
erated on large scales due to non-linear interactions be- 
tween magnetic- and velocity- fields on smaller scales. 
The evolution of magnetic fields arc given by 



5tB = V X (v X B) + 77AB . 



(15) 



We are not interested in the detailed magnetic field evo- 
lution, but rather in the average buildup of magnetic 
energy density taken over an ensemble of cosmic realiza- 
tions, i.e. 



{dtMk) 



(16) 



with Alk the magnetic spectral energy (not to be confused 
with the magnetisation) defined through 



(17) 



where V denotes volume and where we have implicitly 
assumed cosmic homogeneity and isotropy rendering Mk 
a function only of the magnitude of wavevector k. 

Ensemble averaging over the product of Eq. ()15p and 
B, assuming time-independent and uncorrelated fluid 
turbulence, i.e. (wiJ(t)wk' (^0) ^ (^k'k^t't unmodifled by 
the magnetic fields, i.e. neglecting the backreaction mag- 
netic force term in Eq. (1141 Such a calculation has been 
performed first by Ref. [la]. We follow here the anal- 
ysis of Ref. [l^ which leads to the following evolution 
equation for Mk [l7| 



dt 



C / dfcd6' 



kkl v{Lk) 



fiq/k,0)Mk, (18) 



where for simplicity we have dropped the ensemble av- 
erage brackets and discarded terms of higher order in 
wavevector k as well as terms containing magnetic helic- 
ity. In other words, we are only interested in the small 
k tail of non-helical fields. Eq. (|T5|) describes the gen- 
eration of magnetic energy on scale with wave vector q 
by the non-linear interaction between magnetic modes 
with wave vector k and velocity modes with wave vec- 
tors ki = q — k. Here 9 is the angle between q and k, 
v{Lk) are typical velocities on scale Lk — iir/k associated 
with wave vector fc, and C is a crucial constant estimated 
approximately as C ~ (27r)^. Furthermore, the function 
/ is given by 



(1 



X — xcosv 



(1 



2a;cos^) 



(19) 



where the denominator is ki/k. Wc now assume a power 
law for the scale dependence of the fluid velocity, i.e. 



Lk 



(20) 



where Lq = 27r/fco is some reference scale on which 
V = Vg. Here a = 3 corresponds to uncorrelated white 
noise whereas a = —2/3 corresponds to a red Kolmogo- 
roff spectrum. Inserting Eq. (|20p into Eq. ([THj), changing 
to comoving wave vectors k with fcphys = fc/a, where a 
denotes the scale factor, one finds 



kj{a) 



dlnfc A; 



V, 



k 



1/2 



( ^ 1 Mk 



dMg 

d\na 

^^^^^ 

for the evolution of magnetic energy density during a 
radiation dominated Universe. Here fc/(a) is some effec- 
tive ultraviolet cutoff, or damping scale, to be discussed 
below, and Hg is the Hubble constant at the initial mag- 
netogenesis era where we have chosen scale factor a = 1. 
Note that equation Eq. ([^ really holds for some prop- 
erly defined comoving Mk which takes into account the 
redshifting ^ 1/a* of magnetic energy density [isj . 
During matter domination the dependence on the scale 
factor of Eq. (j^TI) is changed from a to a^/^. The function 
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G" in Eq. (1211) contains the 6 integration 



G"(a;) 



1/1 • 3//,\ (1 + a: - a;cos6') 

de sm-' (6*) ^r -^t:^ — ^ 

^ ' {l + x^ ~ 2a;cos6l 5/2-"/2 



(22) 



We are only interested in large-scale modes q'Z k. In this 
case G" (x) is a slowly varying function, with, for example 
varying only from 4/3 to 5/3 when x goes from to 
1. We therefore conservatively drop it. 

Eq. (|18p represents our master equation. It demon- 
strates that magnetic energy may be transferred from 
small scales (large k) to large scales (small q) via non- 
linear processes between fluid flows and magnetic fields, 
and that this process apparently becomes more efficient 
at later epochs (larger a). Before presenting a prediction 
for the evolution of the large-scale tail of magnetic fields, 
we have to emphasize a conceptual difference between 
the findings of this work and all other prior studies of pri- 
mordial magnetic fields. So far, it had been believed that 
the epochs of magnetogenesis and subsequent field decay 
are well separated. In other words, it would be assumed 
that some magnetogenesis scenario during a cosmic phase 
transition would fully account for all field properties, in 
particular, also those on large scales, and that the final 
late time field would be simply given by all magnetic en- 
ergy which hadn't been dissipated yet. This study shows 
that magnetogenesis should be considered as a continous 
process, which in particular holds true for the large-scale 
tail, with magnetic fields generated on small q even long 
after the epoch of the initial magnetogenesis era. 

We now proceed to estimate the growth of large-scale 
magnetic fields. We assume that the kinetic energy den- 
sity is peaked on some scale fc/, which we refer to as the 
integral scale. For k > kj the velocity spectrum is red 
(e.g. Kolmogoroff), and for k < kj given by Eq. pO]). 
Contributions from v modes with fc > fc/ to the integral 
Eq. may be approximatively neglected, as already 
done in Eq. (PT|) . We assume that the magnetic seed 
field only exists on some small scale fc^, with Mq = for 
either q ks and q^ ks- For the moment we make the 
crucial assumption that ks ^ kj and that equipartition 
holds on that scale, i.e. Mk « Ek, with ex = j dkEk, 
where sk is the kinetic energy density. We will discuss 
this assumption in some more detail below. 

We need a model of the evolution of fluid flows in the 
expanding Universe, as this modifles the v{Lk) entering 
Eq. ( 21 1 . Such a model has been presented, in detail, in 
It had been found that the evolution of velocity 
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Ref. 

and magnetic fields in the early Universe is quite com- 
plex, described by an alternation of turbulent phases and 
viscous phases. Here, early on, viscosity is due to neu- 
trinos, and at later stages, due to photons. Nevertheless, 
we focus first on turbulent phases. 



A. Turbulent Phases 

In the absence of further turbulence producing mech- 
anisms, turbulent flows simply decay when the eddy 



turnover rate equals the Hubble rate, i.e. 

yPiki) _ vP{ki)kP 



2-K 



H . 



(23) 



In particular, eddies on kj break up into eddies of larger 
k, with their energy getting transported to the dissipa- 
tive scale fcdiss > kj where it is converted to heat. Here 
the superscript p denotes proper quantities. These are 
related to comoving quantities via La, k^ = k/a, 

as well as = V (radiation dominated) and = v/a 
(matter dominated). From Eq. using Eq. (|20p and 

assuming radiation domination one derives 



fc/(a) = fc?a-2/("+2V 



ith 



2TTHn 



(24) 



(25) 



the initial integral scale. Here we adopted the reference 
scale fco = kj such that ef^ = p-^Vg/2 represents the initial 
turbulent energy density. 

Assuming equipartition Mk{k'^) sa E^^kj) on the initial 
integral scale kj but Mk ~ beyond, we are interested to 
know if the decaying magnetohydrodynamic turbulence 
may establish equipartition on even larger scales, and 
if yes, how large is the largest scale where equipartion 
may be achieved. Using Eq. (|24p . assuming equipartition 
Mk ~ Ek on ki at all times (to be verified selfconsistently 
below) , one may derive for the change in magnetic energy 
density AMq{a) at redshift a 



AMq{a) ~ Cpr 



2TTHn 



4 7 

Ho 



(26) 



where pr is radiation energy density and where we ne- 
glected an unimportant numerical factor. The total 
change of Mg is conservatively determined at Ug where 
the integral scale has decreased to q, i.e. ki{aq) = q. The 
condition that by Uq the change AMq{aq) > M^^i"' = 
Vg{q/kj)°'/q is larger than the equipartition value then 
yields 



C: 



1 

2^ 



(27) 



independent (!) of q, Vg, Hg, and a. As we determined 
C w (27r)2 thus equipartition between magnetic fields and 
velocity fields may be attained on a large range of scales, 
independent of the turbulence spectrum. This justifies 
selfconsistently our above assumption. 

At any late time, the magnetic field spectrum therefore 
has the form 



Mqja) 
Mj{a) 



ki(a) 
Q 

fcj(a) 

fccqui(a) 
ki(a) 







-1 






q 







ki{a) <q< fcdiss (a) 



fcoqui(a) <q< ki{a) 



q < fccqui(a) , 



(28) 
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with the g— independent constant Mj{a) = M^'^^^y This 
is Kolmogoroff (/3 = 2/3) or Ironshnikov-Kraichnan {(3 = 
1/2) on very small scales, following the fluid spectrum 
on scales larger than the integral scale kj (a) and turning 
over to a steeper spectrum for k < fcequi(o)- Here fccqui(a) 
is the smallest wavevector for which equipartition could 
be attained at rcdshift a. It is given by 

A:equi(a) =^ (2^C)^/c/(a). (29) 

The spectral index 4 for fc < fccqui in Eq. pS]) is dic- 
tated by the power of q in Eq. (PT|) . The evolution of the 
spectrum for a reprcscntivc case is displayed in Fig. [TJ 
with the spectrum shown at temperatures T — 100 MeV, 
T = 100 keV, and T = 1 eV, respectively. Note that since 
the integral in the master equation (|2ip is dominated by 
the integral scale, these results are insensitive to the tur- 
bulence power spectrum at scales k <^ kj{a), and thus in 
particular on whether this spectrum is causal or not. 

Eq. (|24[) and Eq. (j26p apply only during radiation dom- 
ination. It is easy to show [l3| that during matter dom- 
ination ki{a) stays constant or at best decreases loga- 
rithmically with scale factor. Using this plus the matter- 
domination modified Eq. ([2T|) . as well as taking into ac- 
count of the proper redshifting of v, it can be shown 
that no further substantial growth of Mg occurs during 
matter-domination. This holds true even when v va, 
as should be the case, yielding a milder redshift depen- 
dence V ^ a~^/^ than the above stated v ^ a~^. Here va 
is the Alfvcn velocity. This implies that fccqui gets " frozen 
in" at the transition of radiation to matter domination 
at its minimum value fc™™;. 

We comment here shortly on one implicit assumption 
made. In particular, we allow Mg only to grow to its 
scale-dependent equipartition value Mq « Eg. Though 
this seem fairly plausible, it is not clear if, in fact, it is 
true. Rather, it is possible that on scales q < kj the mag- 
netic field grows to super-equipartition, i.e. Mq > Eg. 
This may be possible as long as global energy conserva- 
tion is not violated. It is not clear if this could happen, 
and only an improved analysis taking account of potential 
modifications of Eq. (fTSl) due to backreaction of magnetic 
fields on fluid flows, as well as deriving an independent 
equation for (dtEk), could conclusively adress this issue. 
This formidable task is however beyond the scope of this 
paper. We note that our assumption Mq ^ Eg is conser- 
vative. 



B. Viscous Phases 

A realistic calculation has to take into account the vari- 
ous small Reynolds number R < I viscous phases of mag- 
netohydrodynamics in the early Universe. Here viscous 
forces in Eq. (|T4)) are given by [l3| 

f„ = ( ^""'^ t ' (30) 
[ -av /nifp :S> L ' 




10° 10=^ 10* 10" 10" 

k_comoving [Mpc '] 




10° 10^ 10* 10*^ lO' 

k_comoving [Mpc^^] 




10° 10^ 10* 10'^ 10' 

k_comoving [Mpc ^] 



FIG. 1: Evolution of the magnetic energy spectrum (in units 
of the radiation energy density) according to Eq. (|2H) for a 
magnetogenesis temperature T = 100 MeV and for Vg = 1 and 
a — 3. The upper panel shows the assumed initial condition 
for the magnetic- (solid) and kinetic- (dotted) energy spec- 
trum at r = 100 MeV, whereas the middle and lower panels 
show these quantities evolved to T = 100 keV and T = 1 eV, 
respectively. Dashed and solid vertical line represent the inte- 
gral scale at the initial and current temperature, respectively. 
The scale at which the solid line departs from the dotted 
line, is the scale fcoqui from Eq. (|29p which in the present case 
reads fcequi(a) — 0.063fcj(a). This is nicely confirmed by the 
simulation. Note that the semianalytic simulation assumes a 
turbulent regime throughout. 
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where Zmfp is the mean free path of the momentum trans- 
porting particles, early on neutrinos, and later photons, 
and L is the scale under consideration. It can be shown 
that the diffusion constant rj and drag constant a are 
essentially rj ~ /^fp and a ~ 1/^mfp up to numerical, 
and statistical weight factors. Starting at high tem- 
peratures, as the mean free path of neutrinos increases 
with the expansion of the Universe, an initially turbulent 
fluid becomes viscous due to neutrino diffusion. At this 
point all pre-existing fluid motions will be erased up to 
the neutrino diffusion length, i.e. d ~ ^yrft with t cos- 
mic time. However, importantly, magnetic fields are not 
erased, and remain frozen into the plasma. On the other 
hand build-up of large-scale magnetic fields via Eq. (fT8|) 
has stopped entirely as there are no appreciable fluid ve- 
locities. This viscous phase occurs already at fairly high 
temperatures T ~' 20— 80 MeV, depending on the kinetic- 
and magnetic- energies. As the neutrino mean free path 
increases beyond the scale Lj ~ l/kj, the fluid continues 
its viscous period in the neutrino free-streaming regime. 
Form their on, the very small Reynolds number starts 
increasing again. Here secondary fluid velocities 



ki 
pa 



aLr 



(31) 



are generated via the magnetic stresses in Eq. p^ . These 
velocities become larger and larger, and at some point, 
still in the viscous regime, the "instantaneous" dissipa- 
tion scale fcdiss falls below the integral scale kj at the 
beginning of the viscous regime. At this point further 
magnetic- and fluid- energy density is dissipated lead- 
ing to an increase of fc/ « fcdiss- Here fcdiss is given 
again by the condition Eq. (j23p with velocities now as in 
Eq. Somewhat later, typically at T 5 - lOMcV 

shortly before neutrino decoupling, the fluid becomes 
again fully turbulent. This cycle repeats itself at tem- 
peratures T ^ 1 MeV, but now due to the long mean free 
path of photons. 

It is known [l^ that the net effect of theses viscous 
phases is a delay of the growth of the integral scale Li (a) 
during viscous periods themselves. Nevertheless, once 
turbulence has recommenced, the integral scale Li{a) is 
indistinguishable from LY^'°{a), the respective scale in 
a Universe which had always remained turbulent. This 
may be observed in Fig. [5J where results for fc/ (a) of a 
semi-analytic calculation of L/ are shown. The simula- 
tion takes proper account of the correct neutrino- and 
photon- viscosities. It is far from clear if primordial 
viscosity has a similar effect on Lcqui(a) the maximum 
scale where local cquipartition on large scales has been 
attained, in particular, simply the delay of growth of this 
quantity. During large parts of the viscous eras no growth 
of magnetic fleld energy density via Eq. ([T5)) may oc- 
cur. This could yield much weaker Mq. On the other 
hand, towards the end of the viscous periods Mkj en- 
tering Eq. (jl8|) is genuinely larger than when viscosity 
would be small. It turns out that both effects roughly 
cancel each other as can be seen in Fig. [21 The detailed 
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FIG. 2: Evolution of various length scales as a function of 
temperature. All length scales are in units of the approx- 
imate QCD Hubble radius, i.e. if^^ ~ 1 pc. Shown are 
the integral scale Li on which most of the magnetic- and 
kinetic- energy is contained, as well as the cquipartition scale 
icqui, i.e. the largest scale on which equipartition between 
kinetic turbulent energy and magnetic energy could be at- 
tained. Each length scale is shown twice, once for a Universe 
which is always turbulent (dotted lines) and once for a realis- 
tic Universe, where turbulent and viscous MHD are alternat- 
ing due to neutrino- and photon viscosity (solid lines). The 
dashed lines, show the neutrino- and photon- mean free path 
as labeled. The results of this semi-analytic simulation as- 
sume an initial magnetogenesis epoch at T = 100 MeV with 



= 10 ■'pj? and = 10 ^pp, where efj, e^, and are 



initial magnetic-, kinetic-, and total radiation energy density, 
respectively. A white noise a = 3 spectrum has been assumed 
for the slope of the turbulence on large scales Lj. 



simulation shows that as in the case of ki[a), kequiio.) 
during turbulent phases is as if viscous phases had never 
occurred. 



We now briefly comment on the adopted initial con- 
ditions of equipartition on kj, though we do not in- 
tend to go into much detail of the initial magnetoge- 
nesis period itself. Using Eq. ([T8)) one may derive an 
approximate equation dM^Jdt ~ Mkjr with r given 
by « Cv{kj)kj and thus is a fraction of the eddy 
turnover time on the integral scale. We note that since 
v{kj)kj / (2tt) ^ H the requirement to have more than one 
e-fold of amplification of magnetic fields on kj per Hubble 
time during the phase transition reads C^2tt. Since this 
is the case, even initially much weaker than equipartion 
fields are likely to come to equipartion during the phase 
transition. Nevertheless, this does not adress concerns 
about possibly most magnetic energy going into large q 
modes, due to the stronq g"* dependence in Eq. (jl8p . 
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C. Results 

Given all the above, the final present-day spectrum, 
barring further evolution in the non-linear structure for- 
mation regime, is therefore given by Eq. (j28p with fccqui 
and ki replaced by A:™"; and fc™'", their respective values 
at matter-radiation equality. In this subsection we devi- 
ate from our above chosen convention that a = 1 at the 
initial magnetogenesis epoch. Rather, as customary, we 
choose a = 1 at the present epoch. In this case kf™ and 
Mj{aeq) in Eq. (^5)) are given by 



kf' = kiM 
Miiaeq 



2t:H„ f H, 



V 



Q + 4 
2q + 4 



) = Pr^q -TT- 



ki{a, 



(32) 
(33) 



eql 



where Heq and a^q 



5 X 10'^ are Hubble constant 
and scale factor at matter radiation equality, respec- 
tively. We note that for Vg ~ 1, i.e. the speed of light, 
their is equipartition between kinetic turbulent energy e^j- 
and radiation density (as well as subsequently mag- 
netic energy density eg). Such equipartion corresponds 
to Bq ~ 5.7 X 10~^ Gauss in comoving magnetic field 
strength. Making the reasonable assumption that the 
turbulence is white noise (cf. [l9| for a numerical justi- 
fication), i.e. a = 3 yields spectral index n = 3/2 for 
k > fccqui and only for fc < fcoqui the in Ref. @ proposed 
n — 5/2. We will see that this leads to strongly modi- 
fied predictions for the possible large-scale magnetic field 
strength arising due to "causal" processes compared to 
those given in Ref. [I]. The scale L^™- = 27r/fc, 
be obtained via Eq. ([59]) and Eq. 



cqui 



may 



L^-i^30kpc 



1/2 



VlOOMeV 



-3/5 



for 



To derive an upper bound we assume the maximum 
plausible energy in the turbulence, i.e. equipartition 
~ pf. with the radiation energy density. We derive 
the magnetic field strength on scale L from Eq. ([25)1 
and Eq. ((32)) - Eq. (p3)) taking account of the conver- 
sion given below Eq. ([55)1 . On L ~ 100 kpc, the scale 
chosen by Ref. [3 we find B(lOOkpc) < 5 x lO^^"* Gauss 
for the QCD transition (Tg ~ 100 MeV). This value is 
also refiected by Fig. (T] lower panel when using that 
Bk ~ 5.7 X 10"*^ {Alkk/prY/"^ Gauss. For the electroweak 
transition (Tg ~ 100 GeV) one has S(100kpc)<3 x 
10~^° Gauss. This should be compared to the much 
smaller B(lOOkpc) < 10"^^ Gauss given in Ref. Q for the 
electroweak transition. 

However, we stress here again 0, [l^ that more im- 
portant than the magnetic field strength on some large 
scale ~ 100 kpc should be the total magnetic energy on 
all scales. Numerical simulations show [2^ that gravita- 
tionally induced inverse cascades during cluster collapse 
seems to transfer this magnetic power to large scales. 



We thus restate here [1, [T^ the prediction for present 
day magnetic field strength and coherence length as a 
function of magnetogenesis temperature Tg, initial mag- 
netic/turbulent energy density, Tg (in units of radiation 
energy density), and initial helicity hg (as a fraction of 
maximal helicity) 



L.{T) . 12 PC (^) 



1/2 



T 



100 GcV 



-3/5 



Bc{T) ~ 6.0xlO-"G /^V^' ( ^'jl 



0.01 



100 GeV 



-3/5 

(35) 



for non-helical fields and 



Lc(r) ~ 1.9 kpc 



Ta \^/^ f h 



B^{T) ~ 1.6 X lO^^'G 



0.01/ VO.Ol 

-11 



1/3 



100 GeV 



-1/3 



(36) 



rq V'^ ^ hg 



0.01 



0.01 



1/3 



100 GeV 



-1/3 



for helical fields. Here Eq. (|35)) assumes a spectral index 
n = 3/2. 

V. CONCLUSIONS 

In this paper we have considered the large-scale (small 
wavenumber fc) tail of primordial magnetic fields pro- 
duced "causally" (i.e. not during an infiationary epoch) 
in the early Universe. Here we focussed on scales which 
are instantaneously super-Hubble during the initial mag- 
netogenesis epoch, i.e. A 1 pc, as such scales correspond 
to current galactic scales. It has been claimed @ that due 
to causality such fields are highly supressed. i.e. n = 5/2 
in Eq. ([1]) . We have shown here that this may not be the 
case. 

It is customarily assumed that an initial period of mag- 
netogenesis (e.g. during a phase transition) may be well 
separated from the subsequent epoch of freely decaying 
MHD. In other words, field production only happens dur- 
ing the magnetogenesis epoch, whereas subsequently only 
field dissipation takes place. However, this is not the 
case, magnetic field on large scales may be continously 
generated due to nonlinear interactions of small scale 
fluid- and magnetic- modes, even well after the initial 
epoch of magnetogenesis. We have shown that, with the 
expansion of the Universe, this process is rapid enough to 
yield scale-dependent equipartition between fluid veloci- 
ties and magnetic fields on ever and ever larger scales. In 
particular, barring super-equipartion values of the mag- 
netic field, this implies that the kinetic- and magnetic- 
energy tail could have the same spectral index. As white 
noise rt = 3/2 for the fluid is expected [l^, the magnetic 
field tail could also develop n = 3/2. 

When the realistic alternation of turbulent and vis- 
cous MHD (due to photon- and neutrino- fluid viscos- 
ity) in the early Universe is considered, the evolution be- 
comes substantially more complicated. We have followed 
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this with a semi-analytic treatment. During turbulent 
epochs the fluid forces equipartition with magnetic fields 
on small scales. After this, amplification of large-scale 
magnetic fields due to non-linear interactions between 
velocity- and magnetic- modes on small scales becomes 
important. When fluid viscosity becomes considerable, 
all fluid flows up to the neutrino (photon) diffusion scale 
become erased, and the process of amplification of large- 
scale magnetic fields effectively stops. Only when neu- 
trino (photon) interaction with the plasma become very 
weak, new fluid flows are generated by the priorly gener- 
ated and frozen in magnetic fields. Thus the fluid initial 
kinetic energy had been stored in magnetic fleld energy. 
Subsequently further large-scale magnetic fleld amplifl- 
cation takes place. As complicated as is, the net result of 
these processes is a present day large-scale magnetic field 
tail essentially indistinguishable from that if the Universe 
would have stayed turbulent throughout. 

However, we caution that our results are preliminary, 
as our evolution equation may not include all magnetic 
sink terms on large scales. Those were assumed to be- 



come important only at scale-dependent equipartition, 
i.e. when Mq sa Eq. Furthermore, backreaction of mag- 
netic fields on fluid flows was only taken into account in 
a heuristic way. A complete analysis, which is beyond 
the scope of the present paper, should be performed to 
conflrm our results. 

The resulting magnetic fleld spectrum may thus be n = 
3/2 up to some large scale in Eq. (p4| and n = 5/2 
beyond. This would lead to considerably larger magnetic 
field strength on large scales than claimed in Ref. 0. 
If appreciable primordial magnetic fields actually exist, 
however, can be verified only by measuring the magnetic 
field in the intergalactic medium. 
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